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Decays of the Xcj states (J=0, 1, 2) to AAtt + tt~ , including processes with inter- 
mediate £(1385), are studied through the El transition tp' — > 'yxcj using 106 million 
tp' events collected with the BESIII detector at BEPCII. This is the first observation 
of Xcj decays to the final state AA7r + 7r~. The branching ratio of the intermediate 
process Xcj — > S(1385) ± S(1385) T is also measured for the first time, and the results 
agree with the theoretical predictions based on the color-octet effect. 

PACS numbers: 13.25.Gv, 13.30.Eg, 14.20.Pt 
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I. INTRODUCTION 



Decays of P-wave charmonium states, e.g., the Xcj-, cannot be well explained by the 
color-singlet contribution alone, although this works well in explaining the decays of S- 
wave charmonium, e.g., the J/0 and ip'. In calculations of the color-octet contribution, 
Ref. [1] predicted branching ratios of XcJ — >-baryon+anti-baryon in which the Xcj PP 
result is consistent with experimental observation, while the Xcj AA [2] result is not. 
The calculated branching ratios are B(xd -> AA) = (3.91 ± 0.24) x 10~ 5 and B(xc2 -> 
AA) = (3.49 ± 0.20) x 10~ 5 , while the experimental results are (11.8 ± 1.9) x 10" 5 and 
(18.6±2.7) x 10~ 5 , respectively. In addition to AA, reference [1] also calculated the branching 
ratios of Xa -> £(1385)£(1385) and X c2 -> £(1385)£(1385) to be (2.15 ± 0.12) x 10" 5 
and (3.61 ± 0.20) x 10~ 5 , respectively, but there are no previous experimental results on 
these decay channels. Therefore, it is meaningful to test these predictions experimentally. 
In addition, due to the helicity selection rule, the decay of Xco into baryon-antibaryon is 
expected to be suppressed [3]. 

In this paper, we report measurements of Xcj ~^ AA7r + 7r~ (J = 0, 1, 2) (including the 
intermediate £(1385) resonance), XcJ -> £(1385)^^ + 0.0., and Xcj -> £(1385) ± £(1385) T 
through the El transition if}' — > JXcj, where £(1385) ± — > An ± and A — > pn~ . This work 
is based on a 106 million ip' event sample collected with the BESIII detector at the Beijing 
Electron- Positron Collider II (BEPCII) [4]. Continuum data taken at the center of mass 
energy y/s = 3.65 GeV, with an integrated luminosity of 42.9 pb _1 , is used to study non-tp' 
decay background. 
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II. THE BESIII DETECTOR 



BEPCII [5] is a double-ring, multi-bunch e + e~ collider with collision energies ranging from 
2.0 GeV to 4.6 GeV. The BESIII detector [5] is a general-purpose spectrometer with 93% 
coverage of full solid angle. From the interaction point outwards, BESIII is composed of the 
following: a main drift chamber consisting of 43 layers of drift cells with a space resolution of 
about 135 yum and momentum resolution of about 0.5% at 1 GeV/c; a time-of-flight counter, 
which is comprised of two layers of scintillator with time resolution of 80 ps in the barrel 
part and one layer with time resolution of 110 ps in the end-cap part; an electromagnetic 
calorimeter (EMC), which is comprised of 6240 CsI(Tl) crystals, with energy resolution of 
2.5% in the barrel and 5.0% in the end-cap for a 1 GeV photon, and position resolution of 
6 mm in the barrel and 9 mm in the end-cap; a super-conducting solenoid magnet, which 
can provide a 1 T magnetic field parallel to the beam direction; and a muon counter, which 
is made of 1000 m 2 resistive-plate-chambers sandwiched in iron absorbers. 



III. MONTE-CARLO SIMULATION 



For evaluation of the detection efficiency and understanding backgrounds, a Monte-Carlo 
(MC) simulation framework for BESIII was developed. A GEANT4-based MC simulation 
program, BOOST, is designed to simulate the interaction of particles in the spectrometer 
and the responses of the detector. For the generation of charmonium states, e.g., ip', an event 
generator, KKMC [6, 7], is employed, which handles the initial state radiative correction 
and the beam energy spread. For simulation of the resonant decay, BesEvtGen, based on 
EvtGen [8, 9], is used to realize well-measured processes, while LundCharm [8] is used for 
the unknown possible processes. 
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In the MC simulations for the processes presented here, ip' — > 'JXcj is assumed to be a 
pure El transition, and the polar angle, 9, follows a distribution of the form 1 + a cos 2 9 
with a = 1, —1/3, and 1/13 for J = 0, 1 and 2, respectively [10]. Momenta in the decay 
of Xcj -> S(1385) ± S(1385) T , Xcj -> S(1385) ± Avr^(c.c.) and Xcj -> AAtt+tt^ are uniformly 
distributed in phase space. For the decay mode Xcj S(1385) ± S(1385) =F , an extreme 
angular distribution is used to test the phase space assumption and no significant differences 
in efficiencies are observed. This is because, with the current level of statistics, the detection 
efficiencies of the final states are determined mainly by the detection of the El photons, 
and the angular distributions of the hadrons in the subsequent decays are not dependent on 
their MC decay models. 

IV. EVENT SELECTION 

The candidate events for the decay modes ip' — > 7Xcj — > 7AA7r + 7r _ , with A — > pir, were 
chosen with the following selection criteria: 

(1) Charged tracks, i.e., candidates for n ± , p and p, must satisfy |cos 6 1 ! < 0.93, where 9 
is the polar angle with respect to the beam direction. Particle identification is not used. 

(2) The charged tracks not assigned to any A(A) decay candidates must have their point 
of closest approach to the interaction point within 10 cm along the beam direction and 1 cm 
in the perpendicular plane. 

(3) A common vertex constraint is applied to each pair of charged tracks assumed to decay 
from A/A i.e., pn" and p7r + , and the production points of A/A candidates are constrained 
to the interaction point. 

(4) A photon candidate is a shower cluster in the EMC that is not associated with any 
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charged track and has a minimum energy deposit of 25 MeV in the barrel or 50 MeV in the 
end- cap. 

(5) The total momentum of all final particle candidates is constrained to the initial four- 
momentum of the e + e~ system in a kinematic fit. The events with xlc < 80 are retained; 
for an event with more than one photon candidate, only the one with the smallest xic * s 
kept. 

(6) Backgrounds from the decay ip' — > 7i + 7i^J/tp followed by J/tp — > 7AA are rejected 
by requiring the n + 7i~ recoil mass be greater than 3.108 GeV/c 2 or less than 3.088 GeV/c 2 . 
The background from ip' -» £(1385)£(1385), followed by £(1385) TPtx and £° -> 7A, is 
rejected by discarding events with 7A (7A) mass in the range [1.183, 1.202] GeV/c 2 . 

V. SIGNAL ESTIMATION 

The invariant mass distributions of M p7T - and M p7r + are shown in Figs. 1(a) and (b), where 
the signals of A and A are clean. Figure 1(c) shows a scatter plot (M pn - versus M pn +). 
Events where M pn - and M p7T + fall within the box in Fig. 1(c) are used for further analysis. 
The invariant mass distribution of AA7r + 7r - , M AAn + n -, is shown in Fig. 1(d), and the three 
Xcj peaks are clearly observed. 

The invariant masses of Air + and Aix~ (Air~ and Air + ) are displayed in Figs. 2(a) and 
(c) (Figs. 2(b) and (d)), respectively. £(1385) peaks are clearly seen. S 1 * 1 peaks are also 
seen in Figs. 2(c) and (d), around its nominal mass 1.322 GeV/c 2 [2]. The S is a relatively 
long-lived particle, and the selection criteria in this analysis are not optimized for a study 
of the H. Hence, this work does not include study of processes involving S. Events around 
the S 1 * 1 peaks are rejected by requiring M A7T -^ W +^ be less than 1.331 GeV/c 2 or greater than 
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FIG. 1. (a) The invariant mass distribution M pn - for pn~ . (b) The invariant mass 
distribution Mp n + for pn + . (c) The scatter plot of M p7T - versus Mp n +; the box indicates the 
AA signal region used in this analysis, (d) The invariant mass distribution M A ^ 7r + 7r - for 
AA7r + 7r~; the shaded histogram is the background estimated from the inclusive decays of 
the if)' MC sample. 

1.312 GeV/c 2 . 

We divide the remaining Xcj decays into five processes: (1) AAtt + it~' (non-resonant); 
(2) £(1385) + A7T-+c.c.; (3) £(1385)~A7r++c.c.; (4) £(1385)+£(1385)~; and (5) £(1385)~£(1385) + . 
To study the five processes, requirements on M An ~^ w +^ are implemented as shown in Fig. 2. 
The areas between 1.32GeV/c 2 and 1.46GeV/c 2 (two solid arrows) are defined as £(1385) 
signal regions, while the areas smaller than 1.30GeV/c 2 or larger than 1.50GeV/c 2 (two 
dashed arrows) are defined as non-£(1385) regions. 

Due to the broad width and the long tails of the £(1385), the £(1385) and non-£(1385) 
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FIG. 2. The invariant mass distributions of (a) Att + , (b) A7r~, (c) Air~ and (d)A7r + . The 
areas between the two solid arrows are taken as the £(1385) signal regions, while the areas 
outside the two dashed arrows are non-S(1385) regions. The peaks of S 1 * 1 in (c) and (d) will 
be rejected with the requirement M A7r -( A7r +) >1 .331 GeV/c 2 or M A7r -( A7r +) <1.312 GeV/c 2 . 

events feed into the non-E(1385) and £(1385) regions. As a result, the Xcj events that decay 
into the above five processes cannot be completely separated using invariant mass regions 
alone. In this study, we separate the data into five independent categories, with data set 
labels set-j (j — 1, • • • , 5) defined as follows: 



(i) Data set-1: the category to detect the non-resonant process 1. That is, events with 
M\ n +, M^ n ~, M\ v - and Ma 77 + all in non-S(1385) regions. The invariant mass spectrum 
of AA7r + 7r~ is displayed in Fig. 3 (a); 

(ii) Data set-2: the category to detect the single resonant £(1385) + (£(1385)~) process 2. 
That is, events with M A7r +/M A7r - in the £(1385) signal region and with M^-(M A?r +), 



12 

Ma^-, M A7r + in non-£(1385) regions are required. The two types of events in this 
category are combined and displayed in Fig. 3 (b); 

(iii) Data set-3: the category to detect the single resonant £(1385)~(£(1385) + ) process 3. 
Similarly, events with M\ n - /M^ n + in £(1385) signal region and with M^+(M\ n -), 
M& n +, Mj^ n - in non-£(1385) regions are required. The two types of events in this 
category are combined and displayed in Fig. 3 (c); 

(iv) Data set-4: the category to detect process 4. Events with M\ n + , M^ w - in £(1385) 
signal region and M\ n -, M^ w + in non-£(1385) region are selected and displayed in 
Fig. 3 (d); 

(v) Data set-5: the category to detect process 5. Events with M^ n -, M^ w + in £(1385) 
signal region and with M\ w +, M\^~ in non-£(1385) region are selected and displayed 
in Fig. 3 (e). 

The yield in each data set is estimated by a fit to the Xcj peaks, and the yields of each 
process in the full phase space will be disentangled with Eq. (1), as described in Sec. VB. 

The Xcj signal events are clearly observed in each category, as shown in Fig. 3 (a)-(e). 
In the fits of the Xcj m each data set category, a Breit-Wigner function convolved with a 
Gaussian resolution function is used to describe Xcj peaks, while a lst-order polynomial line 
is used to model the background distribution. The Xcj invariant mass parameters are allowed 
to float, while the Xcj widths are fixed to the PDG values [2]. The Gaussian parameters are 
obtained from MC simulation of detector responses. A simultaneous unbinned maximum 
likelihood method is applied, and the fit results are listed in Table I. 
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FIG. 3. The invariant mass distributions of AAtt + 7t~ in the following data samples: (a) 
data set-1, (b) data set-2, (c) data set-3, (d) data set-4, (e) data set-5 and (f) total data 
set. The selections of data set-j (J = 1, • • • 5) are defined in Sec. V. Points with error bars 
are data. The solid curves show the sum of the fitted curves, while the dashed lines are the 
backgrounds. 



A. Background Study 



A 106 million inclusive ^'-decay MC sample is used to investigate possible i/j' decay 
backgrounds. No peaking backgrounds are observed, as shown in Fig. 1(d). Since a large 
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TABLE I. The number of fitted Xcj events in each data set-j (J = 1, • • • , 5) and the total 
data set. nj is the number of fitted Xcj events in data set-j. n to t is that in the total data 
sample. 



Number of events 


XcO 


Xcl 


Xc2 


rti 


10.8±3.8 


12.7±3.9 


36.4±6.4 


U2 


36.4±6.7 


14.7±4.1 


47.6±7.2 


«3 


30.9±6.6 


12.5±4.1 


54.4±7.9 


714 


27.4±5.9 


7.6±3.2 


14.6±4.0 


n 5 


32.8±6.3 


3.6±2.2 


8.7±3.3 


"-tot 


426±23 


105±11 


371±20 



proportion of the Xcj decays are poorly known and their simulations based on the BESIII 
LundCharm model have large uncertainty, we investigate possible underestimated peaking 
backgrounds beneath the Xcj peaks. One major source could be from Xcj AK* + p — > 
pp7r + 7r~K® — > pp27T + 27T~ (c.c); however, the ir + and 7r~ invariant mass distributions of 
candidate events were examined, and no evidence of a K s peak was found. Therefore, 
negligible peaking background is assumed in this study. A study of the continuum data did 
not reveal any non--?// decay backgrounds. 



B. Calculation of Branching Ratios 

To calculate the branching ratios of each mode in Xcj AA7r + 7r _ decay, one has to 
compute the efficiency-corrected number of Xcj decays. The numbers of Xcj events in data 



15 

set-j, which is selected to detect process i, also consists of events from the other processes. 
We describe the number of events of process % in data set-j as 



£11 


£12 


£13 


£l4 


£15 




u\ 




/ \ 

Tlx 


£21 


£22 


£23 


^24 


^25 




N 2 




n 2 


£31 


£32 


£33 


£34 


£35 




N 3 




n 3 


£41 


£42 


£43 


£44 


£45 




N 4 




714 


I £51 


£52 


£53 


^54 


£55 J 




W 




{n 5 j 



(1) 



where Ni is the efficiency-corrected number of events of process i, rij are the numbers of Xcj 
events in the data set-j (as listed in Table I), and denotes the efficiency of process % being 
selected in data set-j, obtained with MC simulation. In practice, the Xcj signals are fitted 
in the five data sets simultaneously, and with the constraint of the three efficiency matrices, 
N\ - N§ are obtained by the fit. Equations (2)-(4) are used to calculate branching ratios 
(£>) of the signal processes, and the results are listed in Table II. The significance of each 
decay mode, which is estimated using Eq. (5), is listed in Table II. Here L m is the likelihood 



of the simultaneous fit, while L 



m(Nj=0) 



is the likelihood of the fit with the assumption that 



Nj is equal to zero. 



B{Xcj — > AA7r + 7r (non — resonant)) 



B(xc.j -> S(1385) +( - ) Att-( + ) + c.c.) 



Nl 



NM^ 1X c.j)B(A^p7T^ 



N- 



2(3) 



NM' 7Xc7)£(£(1385) -+ Att) 



1.0 



B(A -)> p-n-y 



B(xcj -> S(1385) +( ~ ) S(1385) 



-(+h 



N 



4(5) 



7Xcj)^(S(1385) Avr) 2 



1.0 

'B{A ^pnf 



(2) 



(3) 



(4) 
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FIG. 4. The invariant mass distributions of (a) A7r + /A7r~, (b) A7r~/A7r + and (c) 
AA7T + 7r _ . Points with error bars are the data with subtraction of the backgrounds, while 
solid lines are the MC simulation of the signals; The backgrounds subtracted are estimated 
from inclusive MC. The signal components are scaled based on their branching ratios 
measured in this work. The data within 1.312 GeV/c 2 < M^-n^+\ < 1.331 GeV/c 2 are 
removed to reject the S ± candidates. 



Sj — -i 2 x (lnL m — InL 



m(JVj=0) 



(5) 



As shown in Fig. 4, the sum of measured components in the decays of Xcj into the final 
states AA/T + 7r - in MC simulation agrees well with the data. This supports the credibility 
of the decomposition into the different components described above. 



C Xcj -> AA7r + 7r~ (total) 

Based on the selection criteria in Sec. IV, the process Xcj — > AA/T + 7t - (total), including 
the intermediate-resonant processes, is studied. The AA7r + 7r~ invariant mass distributions 
and the fit are displayed in Fig. 3 (f), while the fit results are listed in Table I. According 
to the measured branching ratios of the intermediate resonances in this analysis, signal MC 
samples are generated. This makes the momentum distributions of the final particles in 
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TABLE II. Results of the branching ratios (xlO -5 ) for different decay modes. 'UL' stands 
for the upper limit of the branching ratio at the 90% C.L. ! 5" stands for the statistical 
significance. The first errors are statistical and the second systematic. 



Xcj decay mode 


XcO 


Xcl 


Xc2 


B 


UL 


S 


B 


UL 


S 


B 


UL 


S 


AAtt+tt- (w/o £(1385)) 


28.6±12.6±2.7 


<54 


2.2 


26.2±5.5±3.3 




4.8 


71.8±14.5±8.2 




6.4 


E(1385)+A7r- +c.c. 


34.8±13.2±3.4 


<55 


2.2 




<14 


0.3 


23.6±11.8±2.7 


<42 


1.7 


S(1385)~A7r+ + c.c. 


24.6±12.7±2.4 


<50 


1.6 




<14 


0.0 


37.8±11.8±4.4 


<61 


2.6 


£(1385)+E(1385)- 


16.4±5.7±1.6 




3.1 


4.4±2.5±0.6 


<10 


1.9 


7.9±4.0±0.9 


<17 


2.0 


£(1385)-S(1385)+ 


23.5±6.2±2.3 




4.3 




<5.7 


0.9 




<8.5 


0.0 


AA7r+7r~ (total) 


119.0±6.4±11.4 




> 10 


31.1±3.4±3.9 




> 10 


137.0±7.6±15.7 




> 10 



the MC sample similar to those in experimental data and allows the determination of the 
overall detection efficiency, e t ot: of the sum of all the processes with the same final states 
ifj' — > ^XcJ 7AA7T + 7r _ . The branching ratio of Xcj AA7r + 7r~ (total) is calculated with 
the formula 



B(xcj -»■ AAtt+tt" (total)) 



n to t 



(6) 



VI. SYSTEMATIC UNCERTAINTY 



The systematic uncertainties in this analysis are summarized in Table III. Sources of 
systematic uncertainty include A/A reconstruction, tracking, photon detection, four- 
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TABLE III. Sources of systematic uncertainties. 



Sources 


Relative systematic uncertainty (%) 


A CD 


ACl 


/VCZ 


AA reconstruction 


3.5 


3.5 


3.5 


■k^ tracking (not from AA) 


2.0 


2.0 


2.0 


photon detection 


1.0 


1.0 


1.0 


kinematic fitting 


1.0 


1.0 


1.0 


vetoing background 


4.9 


2.6 


4.4 


fitting method 


4.9 


10.1 


7.9 


the number of tp' 


4.0 


4.0 


4.0 


B(ip' -> 7X c j) 


3.2 


4.3 


4.0 


^(S(1385) ± -> Att*) 


1.7 


1.7 


1.7 


£(£(1385)^ -> Att^) 


1.7 


1.7 


1.7 


total 


9.6 


12.7 


11.4 



momentum constraint kinematic fitting, background rejection, Xcj signal fitting, the number 
of if}' events and branching ratios cited from the PDG [2]. Charged ir tracking and photon 
detection systematic errors are studied following the methods in Refs. [4, 5]. 

For the systematic uncertainty due to A/A reconstruction, J /if; — > AA7r + 7r~ and if>' — > 
7T + 7T~J/tp — > AA7T + 7r _ are used to select a A/A control sample. A/A reconstruction efficiency 
is calculated by taking the ratio of the fitted A/ A yields in the missing mass spectrum before 
and after A/A is found. A/A reconstruction efficiencies consist of tracking efficiency of the 
daughter particles and the vertex-constraint of A/ A. The differences in the efficiencies 
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between experimental data and the MC sample are included in the systematic uncertainties. 

To study the efficiency of the kinematic fitting in the four- momentum constraint, event 
candidates for the three processes ip' — > r yAAn + n^, ijj' — > •"fXcj 73(7r + 7r~) and if)' — > 
J/i/jir + 7r~ — > 3(7r + 7T _ )7r° — > 3(7r + 7r~)27 are used as control samples. The ratio of the event 
rates before and after the kinematic fitting is taken as the efficiency of the kinematic fitting. 
These efficiencies are calculated both in experimental data and in the MC sample, and their 
difference determines the uncertainty of the kinematic fitting. 

For the rejection of the resonances J/V>, and S , different J/ip, S° and S^ 1 

mass region requirements are applied ranging from 3a, 3.5a to 4<r, where a is the detector 
resolution. The largest deviation on the branching ratios is taken as the systematic uncer- 
tainty. The systematic uncertainty of the fitting method is obtained by changing the fitting 
range, the shape of the backgrounds, and changing the detector resolution from the value 
obtained with MC simulation to that obtained by fitting with a free parameter. The relative 
uncertainty of the estimated number of ip' is 4.0% [4]. The uncertainty of the branching 
ratios of intermediate decays are taken from the PDG [2]. The total systematic uncertainty 
is obtained by summing all the individual uncertainties in quadrature. 

VII. RESULTS AND DISCUSSION 

The branching ratios of Xcj decays to £(1385) ± £(1385) =F , S(1385) ± A7r :F + c.c. and 
AA7r + 7r~ (with or without the £(1385) resonance) are measured with 106 million tp' decay 
events collected at BESIII. The results are listed in Table II. The process Xcj AAtt + tt~ is 
observed for the first time. Evidence of Xco ~^ 2(1385) ± S(1385) =F , which strongly violates 
the helicity selection rule, is presented. The branching ratios of Xa,2 — > S(1385) ± S(1385) =F 
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are consistent with the theoretical predictions [1]. 
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